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ABSTRACT 

The Reverse Water Gas Shift process is a candidate technology for water and oxygen production 
on Mars under the In-Situ Propellant Production project. This report focuses on the operation and 
analysis of the Reverse Water Gas Shift (RWGS) process, which has been constructed at Kennedy 
Space Center. A summary of results from the initial operation of the RWGS process along with an 
analysis of these results is included in this report. In addition an evaluation of a material balance 
model developed from the work performed previously under the summer program is included 
along with recommendations for further experimental work. 
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1. INTRODUCTION 


The human exploration of Mars will require the utilization of resources present in the Martian 
environment in order to minimize the payload mass imported from Earth. Reverse Water Gas 
Shift (RWGS), which reacts carbon dioxide and hydrogen to form water and carbon monoxide, 
when coupled with water electrolysis is a candidate technology for oxygen production on Mars. 
The use of the RWGS process for In-Situ Propellant Production, (ISPP) was originally studied by 
Pioneer Astronautics, who determined the RWGS process to be a viable candidate for oxygen 
production. [1 ' 21 

An RWGS prototype system has been constructed and modified at Kennedy Space Center (KSC) 
over the past 18 months. While some experimental runs were made on the RWGS system at the 
end of the year 2000, it was not until July 2001 that successful operation of the process was made 
with no system leaks. Although the operating time has been limited, several conclusions on the 
RWGS operation have been made and problem areas identified as discussed in this report. An 
overview of the process components for RWGS is given below followed by a discussion of the 
results obtained from operation of the system. An evaluation of the validity of a material balance 
simulation model for the RWGS process developed under this research program in 1 999 [3] is then 
presented and followed by a recommendation for further experimental modifications and 
development. 

2. AN OVERVIEW OF THE REVERSE WATER GAS SHIFT PROCESS 

RWGS uses carbon dioxide and hydrogen as reactants to produce oxygen and carbon monoxide 
with a copper on alumina catalyst. Prior to operating the RWGS system the copper catalyst, 
which is produced from copper oxide reduction, must be conditioned to insure no oxygen is 
present in the system. Operating temperatures for processes on Mars are constrained and in 
general designs are not being considered which exceed 500 Celsius. At such temperatures, the 
RWGS reaction equilibrium is lim ited and does not favor the production of water. Thus in order 
to improve the overall conversion a separation step is required to allow the non reacted hydrogen 
and carbon dioxide to be recycled. This is achieved by passing the exit gases to a condenser to 
remove most of the water and then to a hollow fiber polymeric membrane, which preferentially 
permeates hydrogen and carbon dioxide. From the membrane, the reactants can be recycled while 
the byproduct carbon monoxide can be vented. The water produced is stored in a vessel and used 
as feed to an electrolysis unit, which produces oxygen as product and hydrogen, which can also be 
recycled to the reaction process. The RWGS process flow schematic with instrumentation labels 
is presented in Figure 1 below. 

3. RWGS OPERATION AND RESULTS 

The sensors and valves indicated in Figure 1 are referred to for this description of the RWGS 
operation The RWGS process is started up by energizing the heater HR2 and passing a small 
amount of hydrogen (typically 0.1 slpm), through the system from flow controller FC1 or FC2. 
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Figure 1 RWGS Process Flow Schematic 

(Solenoid valves SV2 & SV3) The hydrogen exiting the system is monitored using R9 & RIO 
through the CO vent as the system pressure is regulated between 30 and 60 psia with pressure 
controller PCI. During the start up the recycle compressor is started and the pressure and flow 
are monitored through P8 and R12 respectively. (Measurements P7 and R1 1 are not working) 

It takes approximately 1 !4 hours to heat and stabilize the reactor catalytic bed to 375-380 
Centigrade based on temperature probe T2. A proportional only temperature controller is used to 
aid in achieving this temperature, however it uses T3 as the control variable, which typically 
operates about 35 to 40 degrees cooler that T2 as shown in Figure 2 below. This indirect method 
of controlling the reactor temperature hinders the operation and the time to reach a desired steady 
state due to the manual intervention required. 
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Figure 2 RWGS Reactor Temperatures during July 19 th Operation 

Once temperature T2 exceeds 300 Centigrade, the carbon dioxide flow can be introduced through 
flow controller FC3 or FC4. It was found that waiting until the reaction temperature is closer to 
the desired temperature reduces the start up time since less gas is being heated and the 
endothermic reaction is not started. At this point the water electrolysis can be started if desired 
assuming sufficient water of low conductivity is available in the system. The water is pumped to 
the electrolysis unit using WP1 after passing through an ion exchange bed to purify the water. 

This is a problem area as the exchange resin appears to become saturated after (~1 5 hours of 
operation) and requires new or regenerated resin. 

Another problem with the use of electrolysis in conjunction with the RWGS process is related to 
the control of the system. The largest disturbance is introduced through the water level control 
system. When the electrolysis unit is running, the hydrogen water trap LL1 is continuously 
receiving a water input thus requiring periodic removal. The water is removed by opening 
solenoid valve SV10 whenever the level approaches full. When the valve opens the RWGS sees a 
sudden drop in pressure, which in turn disturbs the flows and temperatures in the system A similar 
disturbance occurs when the RWGS water receiver LL3 is drained after significant water 
production. Another potential problem with electrolysis is the difficulty in maintaining a constant 
source of hydrogen production, since the efficiency of the process is dependent on other process 
variables such as temperature, which are not automatically controlled. 

The RWGS system was operated successfully on 7 different occasions each using a steady state 
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reactor temperature near 379 Celsius. The electrolysis unit was operated the first 3 occasions 
where recycled hydrogen and oxygen were produced. Due to the disturbances discussed above 
however, a sufficient steady state operation could not be achieved and the experimental runs were 
for the next 4 occasions focused on RWGS only with no electrolysis. The 4 th run was successful in 
achieving steady state, however problems with the composition analysis limited the usefulness of 
the data. The best data was obtained during the runs of July 18, 19 and 27 th . 

Table 1 presents a summary of average operating conditions for the runs of July 18, 19 and 27 th . It 
was found from the runs that the overall conversion based on hydrogen, increased with system 
pressure and consequently increased membrane pressure drop. The conversion also increased with 
lower feed rates to the system. It can also be observed from Table 1 that the ratio of product flow 
R9 to H 2 feed decreases as conversion increases since one mole should leave the system for every 
2 moles of reacted feed assuming all of the water is removed at the condenser. 

Steady state operation was closely approached for each of the runs outlined in Table 1 in terms of 
reactor temperature as previously shown in Figure 2 and product flow as shown by Figure 3. It 
can be observed from Figure 4 however, that longer operating times are needed before the recycle 
and membrane feed reach steady state at the higher system pressures. 


Table 1 RWGS July 18, 19 & 27 Operating Summary 


Experimental 

Feed 

mem 

P5 

P8 

Overall 

R9/H2 Feed 

T2 

T3 

T5 

Run 

Flow 

mm 

Port 2 

Port 4 

Conversion 

Ratio 




7/18/01 9:30-10:55 

2.0 

35.94 

52.65 

16.71 

97.37 

1.0146 

379.06 

342.11 

5.67 

7/18/01 1:15-3:45 

mom 

36.71 

52.60 

15.90 

98.66 

1.0175 

378.80 

344.10 

5.61 

7/18/01 4:00-5:25 

mam 

KlETil 

52.76 

18.08 

95.04 

1.0828 

379.82 

342.81 

5.62 

7/19/01 1:15-12:00 


39.42 

59.69 

20.27 

79.46 

1.1508 

378.97 

341.42 

5.88 

7/19/01 12:30-1:55 

■EM 

39.42 

59.69 

20.27 

93.09 

1.0570 

378.97 

341.42 

5.88 

7/19/01 2:15-4:20 


27.11 

44.08 

16.97 

94.62 

1.2870 

378.19 

342.23 

5.75 

7/27/01 12:30-1:35 

2.20 

34.73 

52.63 

ma 

96.37 

1.0559 

378.77 

344.15 

5.46 

7/27/01 3:20-5:05 

2.20 


59.40 

17.34 

99.34 

1.0079 

378.18 

344.92 

5.83 


Table 2 provides a summary of mass balance results obtained from July 1 8 th . A similar table was 
generated for the other 2 days however could not be included due to space restrictions. An overall 
mass balance error between 2 and 7 % was observed over the 3 days of operation. After 
compensating for this error by adding the difference to the flow into the membrane, a component 
mass balance was performed and the resulting error are also presented in the tables. The analysis 
was performed using a gas chromatograph with a thermal conductivity detector for the CO and 
C0 2 and a helium ionization detector for the H 2 . 
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Table 2 July 18 th RWGS Mass Balance Summary 




248 I 297 9155 
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2267 649 


23.17 675 
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2671 3.47 




Pbt2 

Flow 

Rrt 3 
Flow 


220 63.47 1 2861 | 7.92 


235 


245 I 65.67 1 29.82 | 451 


3:00 


3:10 1 63.84 1 31.95 | 421 


320 I 63259 


335 134 1.44 9722 


3:45 130 1.42 9738 




7/1801 


METiTlBli 

415 

425' 

7/1801 

B 

435 

7/1801 

4 

450 

7/1801 

4 

5:00 

7/1801 

IDI 

5:15 

7/1801 

ini 

525 



7142 

2236 

522 

2668 

21.43 

1.89 

7727 

2080 

1.93 

75.76 

2255 

1.69 

75.75 

2257 

L67 


Flow Flow 
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12517 00418 12489 
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1.1771 00415 | 1.1694 


GnpiM F2 CEB 


b2 20864 6785 


CU 3 0036 0029 


Ql4 20564 7.007 


Bitr 0273 -0250 


Rrt 7 fed 
flow flow 



| 2000 

slpm 

1 

00 

girie'mn 

| 2000 

CEBafifited 

1 

00 

giriefrin 


1.960 [ Derail France Brer 
0999 Atkitopcrt2 
1.024 
-0.QZ3 


Rxt7 fed 
flow flow 


27.723 1.018 27.612 


1.1337 Q0416 1.1291 


25.677 1.014 25.486 


1.0500 Q0415 1.0422 



29.612 2000 


12109 I 00818 


27.486 
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36021 
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Q0981 
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1 ZB 
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Figure 3 Product Flow R9 (CO Vent) for July 19 th Operation 



Figure 4 Recycle Flow R12 and Membrane Feed Flow R8 for July 19 th Operation 
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4. RWGS Model Description 

A brief description of the equations used to model the RWGS system is provided in this section. 
The RWGS material balance solution was implemented by programming in MATLAB code. 

A more detailed description of the modeling equations and associated theory can be found in the 
1999 NASA/ASEE Summer Faculty Program Research Report [3] . The chemical equilibrium 
constant K for this reaction is given as a function of temperature by the expression: [4] 

Fro products and reactants and in terms of conversion of a limiting reactant by the expression: 


K = 

5639.5 

e 3 ' 48 r 

, 7 7 49170 

1.077 In T 5.44 x 1 0 4 T + 1 .1 2 5 x 1 0 7 T 2 + ^2 

(1) 

K = 

[CO] [H 2 OJ 

( ©CO + Xeq) ( &H,0 + X e q) 

(2) 

[C0 2 ][h 2 ] 

( 0 C O 2 * Xeq) (1 - X eq ) 


The values of ®co, ®h 2 o, ©C 02 represent the molar ratios of those components to the limiting 
reactant (assumed to be hydrogen) in the inlet stream to the reactor. The limiting reactant is based 
m thermodynamic considerations, the value of K can is related to the concentrations of the on the 
total feed to the reactor, (recycle stream + fresh feed) instead of the fresh feed to the reactor. If 
carbon dioxide were the limiting reactant, equation 2 would need to be modified to reflect this. 
Given a reaction temperature, the value of K can be determined from equation 1 and the 
equilibrium conversion determined by solving equation 2, which is a quadratic in Xe q . The 
equilibrium conversion is the maximum conversion, which can be achieved, in a single pass 
through the reactor. 

The exit gases from the reactor are sent to a condenser where most of the water is removed. For 
modeling purposes, the compositions of the liquid and vapor streams leaving the condenser are 
determined by employing Raoult’s Law to determine the amount of water in the vapor phase and 
Henry’s Law to determine the solubility and hence concentration of the gasses dissolved in the 
condensed phase. 

The gases leaving the condenser are fed into a hollow fiber polymeric membrane which separates 
the components of a gas mixture based on a given components permeability to the polymer. The 
membrane operates in a countercurrent fashion, and is modeled based on a 1 998 publication, ' 
which depicts the membrane as an N stage process as shown in Figure 5. Here, L, k and x J>k are the 
total molar flow and mole fraction of component j in the feed/reject leaving stage k, while V, k and 
yj >k are the total molar flow and mole fraction of j in the permeate leaving stage k. For each stage 
the mass transferred of a component j on a given stage m j k is given by the permeability coefficient 
of that component times the difference in partial pressures in component j across the membrane. 
For N stages, a system of N nonlinear simultaneous equations is yielded for each component j 
requiring an iterative solution. 
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Stage N Stage k+1 Stage k Stage 1 



y j,n y j,N -i y lk+i y j. k y j.k-i y w y j.« 


Figure 5 Flow Diagram for an N Stage Counter-Current Membrane 

After the membrane equations have been solved, the permeate component molar flow rates are 
compared to those assumed initially in the recycle stream and adjusted using the Wegstein 
algorithm until overall convergence is achieved. Due to time constraints, only the model of the 
membrane was evaluated for validity. The errors in the component mass balances resulting from 
the operating data were corrected prior to input into the membrane model. This was done by 
adding or subtracting the required difference needed for mass balance closure to the component 
entering the membrane. The permeability coefficients in the membrane model for each component 
(excluding water) were refined based on a set of operating data on July 1 8 th such that the 
prediction error was zero for the set. The table presented below, shows the model prediction for 
permeate and reject for all the operating runs of July 18, 19 and 27 th . While it can be observed 
that the model predicts well especially in face of the errors in the mass balance, it could be further 
refined with additional data. 

5. CONCLUSIONS 

Successful operation of the RWGS system was achieved during July 2001 confirming RWGS as a 
potential technology for water production on Mars. Various operational and analysis problems 
have been identified to aid in further system development. The model developed in Matlab for the 
RWGS system was validated for the membrane and found to have good predictive capabilities. 
Further testing of the system is needed to better develop the understanding of the RWGS process 
including the examination of the influence of the reactor temperature and the hydrogen to carbon 
dioxide feed ratios. 
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Table 3 Comparison of Experimental Values and Model Predictions for July 18, 19 & 27 



Calculated Actual Absolute Calculated Actual Absolute Calculated Actual Absolute | 

7/18/01 9:30-10:55 

H2 H2 Difference C02 C02 Difference 

GU 

bU 

uinerence 

Feed 

mm m 


0.273 

7.036 

6.786 

0.25 

1.983 

1.96 1 

0.023 

Permeate 


20.564 

0.012 

7.008 

7.007 

0.001 

0.952 


0.072 

Reject 


0.026 

0.013 

0.028 

0.029 

0.001 

1.031 

0.959 

0.072 











7/18/01 1:15-3:45 










Feed 

16.825 

16.861 


7.56 

7.522 

0.038 

2.114 


beeesh 

Permeate 

16.812 

16.812 

0 

7.547 

7.547 

0 

1.126 

mkhm 

BB 

Reject 

0.013 

0.013 

0 

0.013 

HZSBP 

0 

0.988 











— 


7/18/01 4:00-5:25 










Feed 

26.265 

25.333 

0.932 

7.265 

8.126 

0.861 

1.793 


0.071 

Permeate 

26.058 

26.205 

0.147 

7.162 

7.189 

0.027 

0.64 


0.01 

Reject 

0.207 

0.059 

0.148 

0.103 

0.076 

0.027 

1.153 

1.144 

0.009 











7/19/01 1:15-12:00 










Feed 

30.748 

30.555 

0.193 

7.984 

8.07 

0.086 

1.537 

1.644 

0.107 

Permeate 

30.011 

EE 9 

0.469 

7.681 

7.801 

0.12 

0.3831 

■IKW 


Reject 

0.737 

0.267 

0.47 

0.303 

0.184 

0.119 

1.154 


0.158 










— 

7/19/01 12:30-1:55 










Feed 

29.552 


0.124 

8.605 

8.696 

0.091 

2.028 

2.047 

0.019 

Permeate 

29.324 

KsEIifel 

0.139 

8.486 

8.523 

0.037 

0.723 


0.111 

Reject 

0.228 

0.09 

0.138 

0.119 

0.082 

0.037 

1.305 

1.194 











— 

7/27/01 12:30-1:35 










Feed 

25.603 

25.099 

0.504 


7.332 

0.297 

1.986 

1.779 


Permeate 

25.515 

25.563 

0.048 

6.987 

IHsEHSi 

0.027 

0.753 

■IMM 

0.006 

Reject 

0.088 

0.04 

0.048 

0.048 

HB 

0.028 

1.026 

1.032 

0.006 











7/27/01 3:20-5:05 










Feed 

16.507 

17.248 

0.741 

8.921 

IKE2DH 

0.598 

2.47 


IHlufli] 

Permeate 

16.43 

16.5 

0.07 

8.843 

HEE2IDI 

0.066 

1.032 

1.384 

0.352 

Reject 

0.077 

0.007 

0.07 

0.078 

1 0.012 

0.066 

1.438 

1.086 

0.352 
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